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In ultra-peripheral relativistic heavy ion collisions the beam ions scatter at impact parame- 
ters larger than the sum of their radii, so that they interact via long range electromagnetic 
forces. Due to the Lorentz-boost of the beam particles, the exchanged virtual photons have 
high energies and can induce the photoproduction of vector-mesons. We present recent re- 
sults of the STAR experiment at RHIC on p°(770) production in Au-Au ultra-peripheral 
collisions at various energies. STAR has also observed the photoproduction of 7r + 7r~7r + 7r~, 
which is related to the still poorly known excited states of the p°. 

1 Introduction 

The electromagnetic field of a nucleus which is moving at relativistic velocities can be approxi- 
mated by a flux of quasi-real virtual photons using the Weizsackcr- Williams approach [1] . The 
number of photons scales with the atomic charge Z squared, so that fast moving heavy nuclei 
create intense photon fluxes. Relativistic heavy ions may thus be used as photon sources or 
targets. 

In Ultra-Peripheral relativistic heavy ion Collisions (UPCs) the long-range electromagnetic 
interactions are separated from the otherwise indistinguishable hadronic interactions by requir- 
ing impact parameters b larger than the sum of the nuclear radii Ra of the beam ions. Due to 
the large Lorentz-boosts of the beam particles, it is possible to study photonuclear reactions as 
well as photon-photon interactions at high energies in UPCs 0. 

The photoproduction of vector mesons is a typical process in UPCs. A virtual photon, 
radiated by the "emitter" nucleus, fluctuates into a qq pair, which scatters elastically off the 
"target" nucleus and emerges as a real vector meson (cf. Fig.QJi). At high energies the scattering 
can be described in terms of soft Pomeron exchange. The cross section is strongly enhanced 
at low transverse momenta px < 2ft/ Ra of the produced meson, because the qq pair couples 
coherently to the entire nucleus. For these coherent processes the cross section depends on the 
nuclear form factor F(t), where t is the squared four-momentum transfer to the target nucleus. 
For larger pt the qq pairs couple to the individual nuclcons within the target nucleus resulting in 
a smaller cross section which scales approximately with the mass number A modulo corrections 
for the nuclear absorption of the meson. 

Because of the intense photon flux in the case of heavy beam ions, the photoproduction 
of vector mesons may be accompanied by Coulomb excitation of the beam particles. The 
excited ions decay mostly via neutron emission [3] which is a distinctive event signature that 
can be utilized in the trigger decision. In lowest order the vector meson photoproduction 
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Figure 1: Schematic view of the photonuclear production of a p°(770) meson in an ultra- 
peripheral Au-Au collision and its subsequent decay into two charged pions. The meson is 
produced in the fusion processes of a virtual photon 7* and a Pomeron P. a) shows the 
exclusive reaction, b) the one with mutual Coulomb excitation of the beam ions and following 
neutron emission. 



accompanied by mutual nuclear dissociation of the beam ions is a three-photon process. One 
photon produces the vector meson and two additional photons excite the nuclei (see Fig. Uh)- 
In good approximation all three photon exchanges are independent so that the cross section 
can be factorized [3]: 

a V>X nxn = Jd 2 b[l - P had (b)} • P v (b) ■ P xn ,i(b) ■ P xn , 2 {b), 

where -Phad(^) is the probability for hadronic interaction, Pv{b) the probability to produce a 
vector meson V, and P xn ,i{b) the probability that nucleus i emits x neutrons. Compared to 
exclusive photonuclear vector meson production, reactions with mutual Coulomb excitation 
have smaller median impact parameters. 

In this paper we present recent results from the STAR experiment at the Relativistic Heavy 
Ion Collider (RHIC). The Solenoidal Tracker At RHIC (STAR) uses a large cylindrical Time 
Projection Chamber (TPC) [4] with 2 m radius and 4.2 m length, operated in a 0.5 T solenoidal 
magnetic field to reconstruct charged tracks. For tracks with pseudorapidity |?j| < 1.2 and 
transverse momentum pt > 100 MeV/c the tracking efficiency is better than 85 %. The UPC 
trigger is based on two detector systems: The two Zero Degree Calorimeters (ZDCs) [5] which 
sit at ±18 m from the interaction point and measure neutral particles emitted in very forward 
direction. They have an acceptance close to unity for the neutrons originating from nuclear 
dissociation of the beam ions. The second trigger detector system used to select UPC events 
is the Central Trigger Barrel (CTB) [5]. It is an array of 240 plastic scintillator slats that 
surrounds the TPC and provides information about the charged-partiele multiplicity. 

Two basic types of trigger algorithms are used: The "topology" trigger requires a low overall 
charged-partiele multiplicity and subdivides the CTB into four azimuthal quadrants. Events 
with coincident hits in the left and right quadrants are recorded thereby selecting roughly 
back-to-back pion pairs. The top and bottom quadrants are used to veto cosmic rays which 
otherwise could be reconstructed as unlikc-sign particle pairs with zero transverse momentum 
and rapidity. Since there is no requirement on the energy deposit in the ZDCs, the "topology" 
data mainly contain exclusively produced vector mesons. In contrast to this the "minimum 
bias" trigger selects UPC events, where both beam ions dissociated by requiring coincident 
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energy deposits in the ZDCs in addition to a low total charged-particlc multiplicity in the 
CTB. 

In the offline analysis two- and four-prong events arc selected by requiring two and four 
charged tracks, respectively, in the TPC to have zero net charge and to form a common (pri- 
mary) vertex. All tracks are assumed to be pions. In order to suppress backgrounds from 
beam-gas interactions, peripheral hadronic interactions, and pile-up events in addition a low 
overall charged-track multiplicity is required. Backgrounds from pile-up events, beam-gas in- 
teractions, and — in the case of the two-prong sample — cosmic rays are reduced by selecting 
events with their primary vertex in a region close to the interaction diamond. Cosmic ray 
backgrounds in the two-prong "topology" sample are suppressed further by excluding events 
with rapidities y p w 0. Due to the ZDC requirement in the "minimum bias" trigger, the cosmic 
ray background is nearly completely removed. Finally, coherent events are selected by requiring 
small transverse momenta of px < 150 MeV/c for the produced vector mesons. 

2 Coherent Photoproduction of p°(770) 

There arc at least three models that describe the production of p°(770) mesons in ultra- 
peripheral collisions: The model of Klein and Nystrand (KN) [7] uses the Vector Dominance 
Model (VDM) for the virtual photon and a classical mechanical approach for the scatter- 
ing on the target nucleus, based on data from ~/p — > p°(770) p experiments. The Frankfurt, 
Strikman, and Zhalov (FSZ) model [8] employs a generalized VDM to describe the virtual 
photon and a QCD Gribov-Glaubcr approach for the scattering. The model of Goncalves and 
Machado (GM) [5] takes into account nuclear effects and parton saturation phenomena by using 
a QCD color dipole approach. 

The coherent cross section of p°(770) production accompanied by mutual nuclear disso- 
ciation of the beam ions cr c p xnxn is measured using "minimum bias" data. The p° yield is 
estimated by fitting the invariant mass peak of the acceptance corrected m^+^- distribution 
and extrapolating the result from the experimentally accessible rapidity range of \y p \ < 1 to 
the full solid angle using the KN model [7] . In Au- Au collisions at y/s NN = 200 GeV the cross 
section was measured to be cr c ° xnxn = 31.9 ± 1.5 s t a t. ±4.5 sys t.mb; at ^s NN = 130 GeV the value 
is 28.3 ± 2.0 sta t. ± 6.3 syst .mb [TO]. 

Since the efficiency of the "topology" trigger is not well known, the total cross section (Jp ^ ot 
is estimated by applying coherent cross section ratios for different nuclear excitation states, 
which are extracted from the "topology" data, to the cross section values for mutual excitation. 
This way the total coherent p° production cross section is measured to be 530±19 s t a t.±57 sys t.mb 
at = 200 GeV and 460 ± 220 sta t. ± H0 syst .mb at ^/s^ = 130 GeV [10]. In Fig. [§i the 

total cross sections and the cross sections with mutual nuclear dissociation are compared to the 
KN model predictions. 

Figure HJd shows the measured total coherent p° production cross section as a function of 
rapidity and compares to the various model predictions. Due to the limited experimentally 
accessible rapidity range of \y p \ < 1, it is not possible to discriminate the models based on the 
shape of their rapidity distribution. Considering the cross section values, the KN model agrees 
best with the data. 

In p production an interesting interference phenomenon is caused by the fact that the p° 
is produced close (0(1 fm)) to the target nucleus and that the emitter and the target nucleus 
are indistinguishable. Because the impact parameter is larger than the sum of the nuclear radii 
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Figure 2: a) Energy dependence of the total coherent cross section (red) and the one with mutual 
nuclear dissociation (blue) in comparison to the KN model predictions (continuous histogram) . 
b) Comparison of the measured total cross section for coherent p° production with theoretical 
predictions [TU]. The vertical line at each point shows the statistical error. The shaded area 
displays the sum of statistical and systematic errors. The dashed line represents the KN [TJ, 
the dash-dotted line the FSZ [8], and the dotted one the GM model [9|. 



of the projectiles, the system essentially acts like a two-slit interferometer with slit separation 
|6|. Either nucleus A emits a virtual photon which scatters off nucleus B or vice versa. The two 
indistinguishable processes are related by parity transformation and, since the p° has negative 
intrinsic parity, the amplitudes have to be subtracted 

a(p T ,b,y p ) = \A( PTl b,y p )~ AiprA-Vp)^^ 
At mid-rapidity A{px 1 b, y p ) w A{px, b, —y p ) so that the above equation simplifies to 



a(p T ,b,0) = 2\A(p T ,b,0)\' 



1 — cos(pr ■ b) 



The interference is destructive for transverse momenta pt < h/(b)- Figure [3] shows the i(« 
Pj.) distribution, which is roughly exponential at larger t, but has a significant downturn for 
t < 0.0015 (GeV/c) 2 , consistent with the Monte-Carlo simulation that includes the interference 
effect. 

The flight path (3^ct of the produced p° is much smaller than the impact parameter so 
that the p° decays at two well-separated points in space-time. This means that the amplitudes 
overlap and interfere only after the decay and that the interference must involve the tt + it~ final 
state. Interference is only possible, if the final state wave function is entangled, nonlocal, and 
not factorizable into individual wave functions. 

The strength of the interference is extracted from the data by fitting the t distribution with 
the function 



dN_ 
"dT 



— = ae- kt [l+c(R(t)-l)] 
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where k is the slope parameter and c the 
spectral modification parameter that mea- 
sures the interference. A value of c = 
would correspond to no interference, a 
value of c = 1 to the interference pre- 
dicted by the KN model 0E]. The devia- 
tion of the t distribution from the exponen- 
tial shape due to the interference effect is 
parameterized by the function R(t) which 
is determined from the ratio of the simu- 
lated t spectrum with and without interfer- 
ence. The measured spectral modification 
parameter of 87±5 s t a t.±8 sys t.% shows that 
the interference is significant |12| . which 
means that the 7r + 7r~ final state wave func- 
tion retains amplitudes for all possible p° 
decays, long after the decay occurred. The 
system is thus an example of the Einstcin- 
Podolsky-Rosen paradox [TO] with contin- 
uous variables momentum and position. 




'o 0.002 0.004 0.006 0.008 0.01 



t ± [(GeV/c) 2 ] 

Figure 3: Uncorrected t spectrum of p° in the ra- 
pidity range \y p \ < 0.5 for the "minimum bias" 
data [12j . The points represent the data. The 
dashed (filled) histogram is a simulation that in- 
cludes the interference effect, whereas the solid 
histogram is a simulation without interference. 



3 Coherent Photoproduction of 7r + 7r 7r + 7r Final States 

Coherent 7r + 7r~7r + 7r~ production in ultra-peripheral collisions accompanied by mutual nuclear 
dissociation of the beam ions was measured in Au-Au collisions at y/s NN — 200 GeV using 
the "minimum bias" data. The transverse momentum spectrum of the neutral four-prongs 
exhibits an enhancement at low px characteristic for coherent production (cf. Fig. |4^). The 
data show a broad peak in the tt + tt^tt + t:^ invariant mass distribution (see Fig. HJd), similar to 
what was seen in earlier fixed-target photoproduction experiments [13j [14] . This peak is usually 
attributed to the excited p° states p(1450) and p(1700). However, the exact nature of these 
states is still controversial. 

The ir + TT~ir + 7r~ invariant mass distribution was fitted with a 5*- wave Breit-Wigner modified 
by a phenomenological Ross-Stodolsky factor [15] : 

/(m) = A ■ r^Y- 2 m }^\ 2r2 + /BG (m) (1) 

V 777, / (777.Q — m l Y + ITLqTq 

The non-interfering background /bg was parameterized by a second order polynomial which 
was extracted from the invariant mass distribution of +2 or —2 charged four-prongs. Taking 
into account the experimental acceptance, the fit yields a resonance mass of 1540 ± 40 MeV/c 2 
and a width of 570 ± 60 MeV. The Ross-Stodolsky exponent has a value of n = 2.4 ± 0.7, 
however, mass and width depend strongly on the value of n. 

Using the acceptance-corrected ir + Tr~ir + ir~ yield from the above fit and the respective 
p°(770) yield from the 7t + ?t~ invariant mass distribution the cross section ratio C7^ xnxn /<j^ nxn 
is estimated to be 13.4 ± 0.8 %, where again the KN model [7] was used to extrapolate from 
the experimentally accessible rapidity region \y\ < 1 to the full solid angle. Using the measured 
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Figure 4: a) 7r + 7r _ 7r + 7r~ transverse momentum distribution: At low transverse momenta the 
four-prong couples coherently to the entire nucleus leading to a strong enhancement of the cross 
section. The hatched filled histogram shows the expected distribution from simulation. The 
background for the coherent part is estimated from +2 or —2 charged four-prong combinations 
by normalizing their px distribution (gray filled histogram) to that of the neutral four-prongs in 
the region of pr > 250 MeV/c (vertical line) yielding the unfilled histogram, b) Invariant mass 



distribution of coherently produced tt 



The points represent the data, the gray filled 



histogram is the background estimated from charged four-prongs. The thick black line shows 
the fit of the modified S'-wave Brcit-Wigncr of Eq. (p} on top of a second order polynomial 
background (thin black line) taking into account the detector acceptance in the region \y\ < 1 
(rising dotted line). The dashed curve represents the signal curve without background. 



coherent p° production cross section cr c ^n xn the tt^tt tt + tt production cross section is 4.3 ± 
0.3 mb. 

Figure^ shows that tt + tt~tt + tt~ events mainly consist of a low mass tt + it~ pair accompa- 
nied by a p°(770). This motivated the Monte-Carlo decay model p' -> p°(770) Jo (600) which is 
used to estimate the acceptance corrections. As can be seen in Fig. [5^ this model reproduces 
the data well. 

In photoproduction on carbon targets the p' was seen not only in the 7r + 7r~7r + 7r~ decay 
mode, but also in 7r + 7r~ final states |14j . Figure^) shows the high mass region of the measured 
m x + T - spectrum. In order to suppress backgrounds, in particular cosmic rays, tighter cuts are 
applied. The data do not show any significant enhancement around 1540 MeV/c 2 . 



4 Summary 

STAR has measured photonuclcar production of p°(770) in ultra-peripheral rclativistic heavy 
ion collisions. The measured cross sections agree with model predictions. STAR also measured 
for the first time the interference effect in p° production which indicates that the decoherence 
induced by the p° decay is small and that the 7r + 7r~ final state wave function is entangled and 
nonlocal. In addition STAR has observed coherent photoproduction of ir + ir~ ir + 7r~ final states 
in UPCs. The ir + ir~ir + ir~ invariant mass spectrum exhibits a broad peak around 1540 MeV/c 2 
and no corresponding enhancement is seen in the m v + % - distribution. The coherent 7r + 7r~7r + 7r~ 
production cross section is 13.4 ± 0.8 % of that of the /9°(770) meson. 
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Figure 5: a) Invariant Mass distribution of two-pion subsystems: The open circles show the 
measured invariant mass spectrum of the lightest ir + ir~ pair in the event. The filled circles 
represent the invariant mass distribution of the it + tt~ that is recoiling against the lightest 
pair. The spectrum exhibits a clear peak in the p°(770) region. The solid line histograms 
show the prediction from simulation assuming the relative 5-wave decay p' — ► p°(770) /o(600). 
b) High mass region of the m T + T - spectrum with tighter cuts applied in order to suppress 
background: The points represent the data. No significant enhancement is seen in the region 
around 1540 MeV/c 2 where the 7r + 7r~7r + 7r~ invariant mass spectrum exhibits a peak. 
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